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1. Introduction 

The emergence of calculable and detectable vacuum energy as a consequence of imposing 
external boundary conditions in quantum field theory was predicted more than 60 years 
ago by Casimir [TJ. Since that time the vacuum energy of fluctuating quantum fields 
that are subject to boundary conditions was studied in various setups (see, e.g., reviews 
in [21 [3j H]). Usually, the boundary manifold is chosen to be noncompact disconnected 
(e.g., two parallel infinite plates, as generically in pQ) or closed compact (e.g., box or 
sphere); see [21 El H] • 

In the present paper, we shall consider the boundary manifold which is noncompact 
connected and has the form of an infinite tube in three-dimensional space. As has 
been first demonstrated by Aharonov and Bohm [5] in the framework of first-quantized 
theory, the magnetic flux enclosed in such a tube affects the properties of quantum 
matter outside the tube. This effect which is named after them has no analogues in 
classical physics and is characterized by the periodic dependence on the value of the 
flux, vanishing at integer multiples of the London flux value, while being maximal at 
the half of the London flux value. 

In the framework of second-quantized theory, one is interested in the vacuum pola- 
rization effects which are induced outside the tube by the magnetic flux enclosed in the 
tube. In particular, the question is, whether the vacuum energy is induced. If this is 
the case, then this may be denoted as the Casimir-Bohm- Aharonov effect (see also [6]). 

It should be noted that initially [5] the Bohm- Aharonov effect was considered under 
the assumption that the transverse size of the tube is zero, which corresponds to the 
singular magnetic vortex configuration. To take into account the finite transverse size 
of the tube was an important task, since in reality a vortex-forming solenoid is of finite 
widthQ. This task was fulfilled at one time, see [3 [8] . Although, unlike the case of a 
singular vortex, the quantum-mechanical problem in the case of a finite-thickness vortex 
is not exactly solvable, a thorough analysis has been carried out, and, in particular, it has 
been shown that the Bohm- Aharonov effect disappears at a sufficiently large thickness 
of an impenetrable magnetic vortex [7]. 

Returning to quantum field theory, and, appropriately, to the Casimir-Bohm- 
Aharonov effect, we note that up to now this effect was considered for the case of 
a singular magnetic vortex only [61 HI [101 [TTj . Therefore the aim of the present 
paper is to make a first step in the study of the dependence on the thickness of an 
impenetrable magnetic vortex. It should be noted that vacuum polarization effects 
which are induced by magnetic fluxes of finite thickness were considered by different 
authors, see [121 [El ELU EE! HE] . However, these authors are concerned with the case 
when there is no boundary at all and the region of the flux is penetrable for the quantized 
matter fields; therefore, the obtained results have no relation neither to the Casimir, nor 

| Also, a real solenoid is of finite length. However, in the case when the width of the solenoid is much 
smaller than its length and the motion of a quantum-mechanical particle is confined to a plane which 
is transverse to the solenoid, the effects of the width prevail over the effects of the length. 
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to the Bohm-Aharonov effects. In the present paper we shall show that, similarly to the 
Bohm-Aharonov effect, the Casimir-Bohm-Aharonov effect disappears at a sufficiently 
large thickness of the vortex. The second-quantized matter will be represented by the 
charged massive scalar field. 

In the next Section a general definition of the vacuum energy density for the 
quantized scalar field is reviewed and a starting expression for its renormalized value is 
given. In Section 3 this value is calculated numerically in the case of (2 + l)-dimensional 
space-time. Finally, the results are summarized and discussed in Section 4. 



2. Vacuum energy density 

The operator of the quantized charged scalar field is represented in the form 

vf(*°, x) = ^;^= [e-^W) «a + ^VaW b{] , (1) 

A 

where a\ and a\ (b\ and 6a) are the scalar particle (antiparticle) creation and destruction 
operators satisfying commutation relations; wave functions ip\(x) form a complete set 
of solutions to the stationary Klein-Gordon equation 

(-V 2 + m 2 )Mx)=E 2 Mx), (2) 

V is the covariant derivative in an external (background) field and m is the mass of 
the scalar particle; A is the set of parameters (quantum numbers) specifying the state; 
E\ = E-\ > is the energy of the state; symbol ^£ denotes summation over discrete 

A 

and integration (with a certain measure) over continuous values of A. 

We are considering the static background in the form of the cylindrically symmetric 
magnetic vortex of finite thickness, hence the covariant derivative is V = d — ieV with 
the vector potential possessing only one nonvanishing component given by 

V v = $/2tt (3) 

outside the vortex; here <3> is the vortex flux and if is the angle in the polar (r, if) 
coordinates on a plane which is transverse to the vortex. The Dirichlet boundary 
condition on the edge (r = ro) of the vortex is imposed on the scalar field: 

^aU o =0, (4) 

i.e. quantum matter is assumed to be perfectly reflected from the thence impenetrable 
vortex. Provided the orthonormalization condition is satisfied, 

d^x^x = (A|A'>, (5) 

the solution to ([2]) and (HI) in the case of the impenetrable magnetic vortex of thickness 
2r takes form 

^ knkz (x) = (2n)- l e ik * z e in ^ n (kr ) x 

X [^n-e$/27r| (^ ) J| n _ e 0/ 27r | (Ax) - J|„,_ e $/ 2 7r| (&?"o) Y| n _ e $ /2rr| (^ r )] , 
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where z is the coordinate along the vortex, 

(kr )]- 1/2 , (7) 

and < k < oo, — oo < &; 2 < oo, n G Z (Z is the set of integer numbers); J^iu) and Y^u) 
are the Bessel functions of order /i of the first and second kinds. It should be noted that 
the vortex can be obviously generalized to d- dimensional space by adding extra d — 3 
longitudinal coordinates to z; then factor (27r) _1 e lfczZ is changed to (2%)~e lkzZ , where z 
is the (d — 2)-dimensional vector which is orthogonal to the (r, yj)-plane in ^-dimensional 
space. 

In general, the vacuum energy density is determined as the vacuum expectation 
value of the time-time component of the energy-momentum tensor, that is given formally 
by expression 

e = (vac| (doV + d y + d ^d ^ + ) |vac) = ^E x iP* x (x) 1> x (x), (8) 

A 

which is ill-defined, suffering from the ultraviolet divergencies: the momentum integral 
corresponding to the last expression in (jHJ) diverges as p d+1 for p — > oo. The well-defined 
quantity is obtained with the use of regularization and then renormalization procedures 
(see, e.g., [3]). As to regularization, one employs conventionally either heat- kernel or 
zeta-function methods (see, e.g., [2]). As to renormalization, it has been shown [T7] 
that, for a specific configuration of a vortex through the excluded region, it suffices, 
irrespective of the number of spatial dimensions, to perform one subtraction, namely to 
subtract the contribution corresponding to the absence of the vortex. This fact owes 
to the symmetry in the problem, being of rather general nature. It is consistent, for 
instance, with a more recent result obtained in a quite different setup in paper |T8j, 
where the Casimir energy per unit length for n non-overlapping parallel cylinders of 
infinite length in three-dimensional space is shown to be directly related (without the 
need of an extra subtraction or an extra counter-term) to the Casimir energy for n 
non-overlapping discs in two-dimensional space. 

Thus, the renormalized vacuum energy density in the case of the finite-thickness 
vortex takes form 



oo 

e ren = {2-Kf- d J d d ~ 2 k z J dkk(k z 2 + k 2 + m 2 ) 1/2 [S(kr, kr ) - S(kr, £x )|$=o] 



(9) 



where, in view of (jSJ), 
S(kr,kr ) = J2P 2 n(kr )x 



(10) 



X [Y\n-e$/2Tr\(kr )J\ n -e®/2n\{kr) - J r | n -e*/27r|(^o)51n-e*/27r[(^)] • 

Owing to the infinite range of summation, the last expression is periodic in flux $ with 
period equal to 27re _1 , i.e. it depends on quantity 

e$ ITe$T . . 
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where [[«]] is the integer part of quantity u (i.e. the integer which is less than or equal 
to u). 

Let us rewrite ffTUl) in the form 



S(kr, kr ) = So(kr) + S^kr, kr ), (12) 

where So(kr) corresponds to the appropriate series in the case of the vacuum polarization 
by a singular magnetic vortex [91 [10J E] : 

So(kr) = [Ju+fW + Jn+i-F(kr)] = dr [J F (r) J_ 1+j? (t) + J. f {t) J X - F {t)\ , (13) 

n=0 ^ 



and Si(kr, kr ) is a correction term due to the finite thickness of a vortex: 

oo 

S 1 (kr,kr ) = 2j2 



T n » v N J„+F(fcr"o)^n+F(A;r) - y n+ F(fcr )J n+ F(A;r) 



n=0 

Jn+l-f(^o)K+l-^(^) 



J n+ i_F(fcr )y ra+ i_ F (fcr) - F w+ i_ F (fcr ) J n+1 _ F (kr) 
J 2 n+l - F (kr ) + YZ +1 _ F (kr ) 



(14) 



n=0 



T 2 (i \ Jn+F(k r ) + Yn+F(k r ) T 2 (, \ ^n+l-F (^ r ) + ^n+l-F (^ r ) 



In the absence of the magnetic flux in the tube we have 
S(kr, /cr )|$ =0 = S + Si(kr, kr ), 

where 

oo 

S = J*(kr) + 2Y,Jn(kr) = l, 

n=l 

and a correction term due to the finite thickness of an empty tube: 
Si(kr, kr ) = 2 



(15) 
(16) 



j n ver /, \ Jo{kr )Y Q (kr) - Y (kr )J (kr / 
Jo{kr )Y (kr) — — ; : T , 0/ , — - \ 



+ 2j2j n (kr )Y n (kr 



n=l 



J*(kr )+Y 2 (kr ) 

J n (kr )Y n (kr) - Y n (kr Q )J n (kr) 
J*(kr ) + Yi(kr ) 



(17) 



J 2 (fcr ) f^r^£\ + 2 £ J^M- ^ + * W 



J 2 (A;ro) + F 2 (A;r ) ' " ^ "" v " ' u ' J*(*r ) + Y*(kr )_ 

Thus, vacuum energy density flH]) depends on F (jTTj) . i.e. it is periodic in flux $ 
with a period equal to 27re~ 1 . Moreover, relation fl9]) is symmetric under substitution 
F — > 1 — F, vanishing at F — > (F — y 1) and, perhaps, attaining its maximal value at 
F = 1/20. Relations <$3§ and (JUJ are simplified at F = 1/2: 

2fcr 



^(/cr)!^^-! = - 

7T 



sin r, 



T 



§ At least, this is certainly true in the case of the singular vortex both for the Bohm-Aharonov [5] and 
the Casimir-Bohm-Aharonov [§J M US EE] effects. 
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Si (At, kr )\$ =ne -i = 

^ + i(fcr )[r 2 + ,(fcr)-J^( fcr )j -2J n ^(kr )Y nH (kr )J n ^(kr)Y n+ ,(kr) (19) 



Since it is hardly possible to evaluate sums in (JHj) and ifPTj) analytically, our further 
analysis will employ numerical calculation. In the following we restrict ourselves to the 
case of F = 1/2 and d = 2, when the expression for the vacuum energy density takes 
form 

oo 

Sren = — [ dk k (k 2 + TO 2 ) 1/2 GUtT, kr ) , (20) 
27T J 



where 



G(kr, kr ) = S(kr, fcr )|$ =7re -i - S(Ar, Ax )|*=o- (21) 



3. Numerical evaluation of the vacuum energy density 

We rewrite (I20p in the dimensionless form 



3^ 

" 2tt 



oo 

| dzz^z*+(^) 2 G{zM, (22) 







where A = r /r, A G [0, 1]. Let us point out some analytical properties of the integrand 
function in (|22|) : it vanishes at the edge of the vortex 

lim G(z, Xz) = 0; (23) 

A->1 

at large distances from the vortex the case of a singular vortex is recovered 

lim G(z, Xz) = S (z)\^ =ne -i - S ; (24) 

at small values of z one gets 

G(^,Az)|^o = -[ln(A)/ln(A^)] 2 . (25) 

Numerical analysis indicates that in the calculation of function G(z, Xz) one can 
use series in (|17|) and f[T9l) with finite limits, namely for calculation G(z, Xz) at point 
z = z' it is enough to cut off the summation limits by n = [[z 1 + 30]]. In this case the 
relative error is 
G{z, Az)| ne ( 0| [[ z+ 3o]]) - G(z, Xz) 
G(z,Xz) 

It can be shown that the envelope of G(z, Xz) is exponentially decreasing function at 
large z, see Fig.l. So, for the finite-thickness magnetic vortex we can compute values of 
dimensionless quantity r 3 e ren (T2"2"|) for different (not very small) values of A. To do this, 
we have to be able to perform integration in ( 122]) with high precision. We make it in a 
following way. 



< 5(A), 5(A) < 10~ 17 , A G [1/10,9/10]. (26) 
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Figure 1. Behavior of G(z, Xz) at different values of A. 



As one can see from Fig. 2, the function G(z, Xz) is negative from z = to the 
first function root at z = Z\ {z\ ^ 0). So, the appropriate integral in f )22|) is negative. 
The subsequent roots are denoted by Z2, z^, etc. Because of decreasing character of 
the envelope function the integral from z\ to z 3 will be positive. It is useful to define 
a period of function G(z, Xz) as an interval between two next to neighboring roots, i.e. 
from Z\ to zz, from z% to Z5, and so on. Then the full integral in ( 1221) will be a sum 
of the negative integral from z = to z = z\ and a multitude of positive integrals 
over subsequent periods. In the case of sufficiently small transverse size of the tube 
(mro < 0.1) the integrals over some finite number of first periods may be negative but 
thereupon they become and remain positive also. 

For small z (z < 20) we make a direct integration of function G(z, Xz) over periods 
using 25 digits of precision in internal computations. 

For large z we make integration for each period separately. To do it we create a 
table of values of function G(z, Xz) for a separate period and replace this function by 
more simple function in the form 

e -bz j{ f z 2\ 

G int (z,Xz) = 2 , sm(kz + j\nz + (f) ), (27) 

Z JDq yZ J 



G(z,Xz) 




Figure 2. The location of roots of G(z, Xz) at A = 0.7. 
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1.5-10 



1-10 



5-10 




Figure 3. r 3 e ren at mr = 3/2 as a function of x. 



< 10" 



(2* 



where sine function ensures that roots of G int (z,\z) coincide with roots of G(z, Az); 
A q (y) and B q {y) are g-degree polynomials, q can be 3,4 or 5; all unknown parameters 
can be found from an interpolation procedure. We allow a relative error of interpolation 
to be 

G int {z,Xz) - G(z,Xz 
G(z, \z) 

for each period. The function Gi n t(z,Xz) can be immediately integrated with the 
required accuracy. In this way we made integration up to z ~ 100/A with absolute 
accuracy up to 1CT 17 . 

With the help of the above procedure we obtain a table of contributions from 
integration over each period, extrapolate this table to infinity, and after that we find 
the full integral in ( )22l) as a sum of the negative integral over first period(s), a multitude 
of positive integrals over periods up to z ~ 100/A and an interpolation term. The 
absolute accuracy of the obtained result is 10~ 13 . It should be noted that nearly 99 % of 
the integral value in (1221) is obtained by direct calculation and only nearly one percent 
is the contribution from the interpolation. 

Dimensionless quantity r 3 e ren ( 1221) is a function of two dimensionless parameters, 
mr and mr. Using the above described procedure, we calculate r 3 e 



3r ~ ren at several values 



1.5-10 




5-10 



0.5 1.0 1.5 2.0 2.5 3.0 



Figure 4. r 3 e ren at mr = 1 as a function of x. 
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Figure 5. r 3 e ren at mr = 1/2 as a function of x. 

oimr as a function of dimensionless distance from the edge of the vortex, x = m(r — r ), 
in the range < x < 3mro- Further increase of the distance from the vortex results in 
a significant increment of computational time, because there the envelope of G(z, Xz) 
fails to be a sufficiently decreasing function as it is at smaller distances. The results of 
our numerical calculations are presented in Figs. 3-7, where r 3 e ren is along the ordinate 
axis and x is along the abscissa axis; solid lines are interpolating the dots that have 
been calculated. 

The typical behaviour of r 3 e ren is clearly illustrated in cases mr > 1 by Fig. 3 and 
Fig.4 The vacuum energy density is zero at the edge of the vortex (at x — 0), starts 
increasing by some power law x a with a > 1, reaches maximum at x ~ 1, and decreases 
at larger distances to zero (probably exponentially as e~ x ). However, as mr decreases, 
the available range of x is shrunk due to above mentioned restriction x < 3mro- In the 
case of mr = 1/2 a maximum at x ~ 1 is clearly seen (Fig.5), and a following decrease 
to zero may be anticipated. In the cases of mro = 10 _1 (Fig.6) and mro = 10 -4 (Fig. 7) 
one may suppose that there will be a maximum at x ~ 1 and a following decrease to 
zero. But, one can be sure for certain from Figs. 3-7 that the vacuum energy density 
decreases to zero as x a with a > 1 at x — > 0. 




Figure 6. r z e ren at mro = 10 1 as a function of x. 
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5-10~ 5 1 ■ 1 4 1.5-1CT 4 2-1CT 4 2.5-1CT 4 3-10 



Figure 7. r 3 e ren at mro = 10 4 as a function of x. 



As to values of the vacuum energy density, they are rapidly decreasing as parameter 
mro increases and becomes more than unity. Namely, the maximal values of r 3 e ren are: 
3.3 -10~ 7 at mr = 1/2 (Fig.5), 2.1 ■ 10~ 9 at mr = 1 (Fig.4) and 2.2 -lO^ 11 at mr = 3/2 
(Fig. 3). These should be compared with much larger values which are already attained 
below maxima in Fig.6 and Fig. 7: 3.5 • 10~ 5 at mro — lO^ 1 and 7 • 10~ 5 at mro — 10~ 4 . 
It should be noted that, in the case of the singular vortex (mro — 0), the maximal value 
of r 3 e ren is (127T 2 ) -1 fh 8.5 • 10~ 3 [9]; the appropriate plot of r 3 e ren as a function of mr 
is taken from [TT] and is presented in Fig. 8. Thus, one may suppose that, in the case 
of the vortex with thickness in the range < mr < 10~ 4 , the maximal value of r 3 e ren 
will be somewhere in the range 10~ 4 -j- 10~ 3 . For more clarity, the results of Fig.6 and 
Fig. 7 are plotted as functions of variable r/r = A -1 in Fig. 9. Note that, as mr falls by 
three orders from 10 _1 to 10 -4 , quantity r 3 e ren changes by factor 2 only. This should 
be compared with r 3 e ren at m = 0, which is plotted as a function of r/ro in Fig. 10; the 
latter plot coincides actually with that corresponding to mr = 10 -4 in Fig.9. It should 
be noted also that at sufficiently small distances from the vortex edge (at r — ro m -1 ) 
the behaviour of r 3 e ren coincides with that in the m = case. 




Figure 8. r z e ren at r$ = as a function of mr. 
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4. Discussion of the results 

We have studied the influence of finite thickness of the impenetrable magnetic vortex 
on the vacuum polarization of the quantized charged massive scalar field. Since units 
c = h = 1 are used, the London flux value is 2ne~ 1 , and we show that induced vacuum 
energy density ([9]) is periodic in the value of vortex flux <3>, vanishing at integer multiples 
of the London flux value (at $ = 27me _1 ) and being presumably maximal at half of the 
London flux value (at $ = n(2n + l)e _1 ). If the vortex thickness decreases, ro — >■ 0, or 
a distance from the vortex increases, r — r — > oo, then the contribution of Si(kr, Zero) 
( |T4|) and Si (At, fcro) ( |T71) to (El) tends smoothly to zero, and the vacuum energy density 
converges with that induced by the singular magnetic vortex. 

Our numerical analysis of the vortex thickness effects has been carried out for the 
case of the vortex flux equal to half of the London flux value; the quantized scalar field 
is confined to a plane which is orthogonal to the vortex. As follows from this analysis, 
the vacuum polarization actually disappears, when the transverse size of the vortex (ro) 
exceeds the Compton wavelength of the scalar particle (m -1 ): the maximal value of the 
induced vacuum energy density falls by two orders from 2.6 ■ lCT 10 m 3 to 1.4 • 10 _12 m 3 
as mro increases from 1 to 3/2. 

This result should be compared with the influence of the vortex thickness on 
the conventional Bohm-Aharonov effect. In the framework of first-quantized theory, 
one considers elastic scattering of a quantum-mechanical charged particle on the 
impenetrable magnetic vortex of thickness 2tq. The incident wave is characterized 
by momentum p, so the dimensionless parameter of the problem is pr . In the long- 
wavelength limit, pr Q — > 0, scattering converges with scattering on the singular magnetic 
vortex [8]. Since the short- wavelength limit, pro — > oo, corresponds to the case 
when quasi-classical approximation is applicable, one would anticipate that the purely 
quantum effect, as is the Bohm-Aharonov one, disappears in this limit. As has been 
shown in [TJ, this anticipation is indeed confirmed, and scattering in the pr — > oo limit 
converges with scattering of a classical point particle on the impenetrable tube, being 
independent of the enclosed magnetic flux. 




f-» 



1.5 



2.0 



2.5 



3.0 



3.5 



4.0 



Figure 9. r 3 e ren at the smallest values of mr as a function of r/r . 
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In the framework of second-quantized theory, one considers the vacuum polarization 
in the background of the impenetrable magnetic vortex. The appropriate dimensionless 
parameter is mro, and, as we have shown in the present paper, the Casimir-Bohm- 
Aharonov effect disappears in the mr — > oo limit, becoming actually negligible at 
mro > 3/2. 

In the case of the singular magnetic vortex, the induced vacuum energy density 
diverges at the location of the vortex j6j [101 [11] . As it has been shown in the present 
paper, this divergence is unphysical, disappearing when thickness of the impenetrable 
magnetic vortex is taken into account: under the Dirichlet condition for the quantized 
field the induced vacuum energy density is vanishing as (r — ro) a with a > 1 at 
the edge of the vortex. Therefore, the vacuum energy which is induced on the whole 
transverse plane, 



is finite, contrary to the case of the singular vortex when it is infinite. Although we are 
unaware of the value of E ren , the maximal value of e ren is estimated to be somewhat of 
the order of 10~ 3 m 3 if mr < 10~ 4 . 

A brief discussion of polarization of the vacuum of the quantized massless scalar 
field is in order. In this case, the induced vacuum energy density is zero at the edge of 
the vortex, starts increasing as (r — r ) a with a > 1 (see Fig. 10), reaches its maximum 
and then decreases with asymptotics (127r 2 r 3 ) -1 P, [TO] [TT]. Induced vacuum energy 
(129]) is finite in this case also. 

The finite-thickness vortex can be formed as a topological defect appearing after 
a phase transition with spontaneous breakdown of the gauge symmetry [19]. Such a 
structure under the name of a cosmic string [20] [21] is currently discussed in various 
contexts in cosmology and astrophysics, see, e.g., j22j [23]. The cosmic string is 
characterized by flux 27re^ 1 , where en is the coupling constant of the Higgs scalar field 
to the string-forming gauge field; the transverse size of the string is of the order of 



CO 




(29) 
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Figure 10 



• r 3 s ren at m = as a function of r/rp. 
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correlation length m^ 1 , where mn is the mass of the Higgs scalar field. Then, as it 
follows from our consideration in the present paper, the cosmic string can polarize the 
vacuum of quantum matter only in the case when the mass of the matter field is much 
less than that of the Higgs field, m <C m-R- For instance, the cosmic string which is 
formed at the grand unification scale can polarize the vacuum of the electroweak theory, 
whereas the would-be cosmic string corresponding to the electroweak symmetry breaking 
has no impact on the vacuum of quantum matter at the grand unification scale. 
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